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LEAF SEQUENCING METHOD AND SYSTEM 



CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] Not applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

[0002] The United States Government has certain rights to this invention pursuant to NIH 
Grant/Contract No. LM06659-03. 

FIELD OF THE INVENTION 
[0003] The invention relates to a radiation emitting device, and more particularly, to a 
system and method for efficiently and more safely delivering radiation treatment to a patient. 

BACKGROUND 

[0004] Radiation emitting devices are generally known and can be used for radiation therapy 
for the treatment of patients. A radiation therapy device generally includes a gantry which is. 
swiveled around a horizontal axis of rotation in the course of a radiation therapy treatment. A 
linear accelerator is located in the gantry for generating a high energy radiation beam for the 
therapy. This high energy radiation beam can be an electron beam or X-ray beam. During 
treatment, the radiation beam is directed towards a patient lying in the isocenter of the gantry 
rotation. 
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[0005) To control the radiation emitted toward a given object, a beam shielding device, such 
as a plate arrangement or a collimator, is typically provided in the trajectory of the radiation 
beam between the radiation source and the patient. A collimator is a computer-controlled 
mechanical beam shielding device which generally includes multiple leaves, for example, a 
plurality of relatively thin plates or rods, typically arranged as opposing leaf pairs. The plates are 
formed from a relatively dense and radiation impervious material and are generally 
independently positionable to size and shape of the radiation beam. These leaves move across the 
tissue being radiated, thus blocking out some areas and filtering others to vary the beam intensity 
and precisely distribute the radiation dosage. 

[00061 A multileaf collimator (MLC) is an example of a multileaf beam shielding device that 
can accurately and efficiently adjust the size and shape of the radiation beam. The size and shape 
of a radiation beam is designed during the treatment plaiming process. This is useful for both 
intensity modulated radiation treatment (IMRT) and three-dimensional conformal radiation 
therapy (3D CRT). 

[0007] Traditional radiotherapy utilizes uniform beams of radiation, producing a uniform 
distribution of dose throughout the irradiated volume, which includes the target volume. This 
ensures the target is adequately covered, but does nothing to avoid often critical surrounding 
structures. With IMRT, the beams of radiation are made to be intentionally non-uniform. In this 
way the dose distribution can be carefully shaped to minimize radiation to surrounding 
structures. 

[0008] Measurement unit efficiency is a commonly used measure of beam efficiency. The 
measurement unit (MU) efficiency is defined as the efficiency with which the incident radiation 
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results dose being in absorbed in the target region of a patient. The consequence of low MU 
efficiency is an increase in leakage radiation that reaches the surrounding tissue of the patient. 
[0009] There are several critical components of a successful DVIRT program. The first is a 
process referred to as "inverse planning." Inverse planning utilizes a mathematical algorithm to 
optimize the intensity of the various beams. This optimization is highly computer intensive. 
[00010] The second critical process is to convert the intensity distributions obtained, often 
referred to cumulatively as a fluence map, into a series of MLC leaf movements. This is referred 
to as "leaf sequencing." Many device-specific factors must be accounted for in this process. 
These factors include radiation leakage through and between the leaves, leaf speed, dose rate, 
and the "tongue-and-groove" effect. 

[00011] IMRT can be performed either while the beam is on, which is referred to as dynamic 
MLC or DMLC delivery, or by turning the beam off while the leaves move to their next position, 
which is referred to as segmented MLC or SMLC delivery. The beam shielding device defines a 
field on the object to which a prescribed amount of radiation is to be delivered. The usual 
treatment field shape results in a three-dimensional treatment volume which includes segments 
of normal tissue, thereby limiting the dose that can be given to the target, such as a timior. The 
dose delivered to the tumor can be increased and the treatment time decreased as the amount of 
dose delivered to the normal surrounding tissue is decreased. Although current leaf sequencing 
algorithms have reduced the radiation level reaching surrounding normal tissue somewhat as 
compared to traditional uniform beams of radiation, these leaf sequences have not provided 
optimal MU efficiency. 
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SUMMARY 



[0001 2] A method of delivering radiation treatment using multi-leaf coUimation includes the 
step of (a) providing a radiation fluence map which supplies a desired radiation intensity profile. 
The fluence map is (b) converted into a preliminary leaf sequence, wherein the preliminary leaf 
sequence minimizes machine on-time and is generated without any leaf movement constraints. In 
step (c), a leaf movement constraint is imposed on the preliminary leaf sequence. In the next step 
(d), at least one constraint elimination algorithm is then applied, the algorithm adjusting the 
preliminary leaf sequence to minimize violations of the constraint while providing the desired 
fluence map and minimized radiation on-time. The method can be applied to SMLC and DLMC 
systems, and can include adjustment for the tongue-and-groove effect. 

[00013] The leaf movement constraint can comprise a minimum separation distance between 
adjacent leaves in the leaf pair. In this embodiment, the applying step (d) can comprise (e) 
modifying at least one leaf pair in the preliminary leaf sequence to form a modified leaf 
sequence, wherein the modifying step comprises identifying and adjusting positions of leaves in 
the preliminary leaf sequence which violate the minimum separation distance to provide at least 
the minimum separation distance, (f) modifying at least one leaf pair in the modified leaf 
sequence to produce a further modified leaf sequence, the further modified leaf sequence 
providing the intensity profile, (g) examining the further modified leaf sequence for violations of 
the minimum separation distance, and (h) iteratively repeating steps (e) and (f) if at least one 
violation of the minimum separation distance is identified in step (g) using the further modified 
leaf sequence as the preliminary leaf sequence to generate a corrected leaf sequence. In this 
embodiment of the invention, the method can further comprise the step of reducing a tongue- 
and-groove underdose. The step of reducing the tongue-and-groove underdose can comprise 
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applying a tongue-and-groove constraint to the corrected leaf sequence, The step of applying a 
tongue-and-groove constraint can comprise the steps of (i) modifying at least one leaf pair in the 
corrected leaf sequence to form a modified corrected leaf sequence, wherein the modifying step 
comprises identifying and adjusting positions of leaves in the prehminary leaf sequence which 
violate a tongue-and-groove constraint, (j) modifying at least one leaf pair in the modified 
corrected leaf sequence to produce a further modified leaf sequence, the further modified leaf 
sequence providing the fluence map, (k) examining the further modified corrected leaf sequence 
for violations of tongue-and-groove constraint, and (1) iteratively repeating steps (i) and (j) if at 
least one violation of the tongue and groove constraint is identified in step (k) using the further 
modified leaf sequence as the preliminary leaf sequence. 

[00014] In the embodiment where the multi-leaf collimation is dynamic multi-leaf 
coUimation, the constraint can comprise a leaf interdigitation constraint. In this embodiment, the 
applying step (d) can comprise (e) modifying at least one leaf pair in the preliminary leaf 
sequence to form a modified leaf sequence, wherein the modifying step comprises identifying 
and adjusting positions of leaves in the preliminary leaf sequence which violate the 
interdigitation constraint, (f) modifying at least one leaf pair in the modified leaf sequence to 
produce a further modified leaf sequence, the further modified leaf sequence providing the 
fluence map, (g) examining the further modified leaf sequence for violations of the 
interdigitation constraint, and (h) iteratively repeating steps (e) and (f) if at least one violation of 
the interdigitation constraint is identified in step (g) using the further modified leaf sequence as 
the preliminary leaf sequence. 

[0001 5] A method of reducing tongue-and-groove underdose during radiation treatment using 
multi-leaf collimation can comprise the steps of (a) providing a radiation fluence map, (b) 
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converting the fluence map into a preliminary leaf sequence, wherein the preliminary leaf 
sequence minimizes a minimum on-time and is generated without any leaf movement 
constraints, (c) modifying at least one leaf pair in the preliminary leaf sequence to form a 
modified leaf sequence, wherein the modifying step comprises identifying and adjusting 
positions of leaves in the preliminary leaf sequence which violate a tongue-and-groove 
constraint. In the next step (d), at least one leaf pair in the modified leaf sequence is modified to 
produce a further modified leaf sequence, the further modified leaf sequence providing the 
desired intensity profile, (e) examining the further modified leaf sequence for violations of the 
tongue-and-groove constraint, and in step (f) iteratively repeating steps (c) and (d) if at least one 
violation of the tongue-and-groove constraint is identified in step (e) using the further modified 
leaf sequence as the preliminary leaf sequence. 

[00016] A system for delivering radiation treatment using multi-leaf coUimation comprises a 
radiation source for generating a radiation beam, a multi-leaf collimator having a plurality of 
leafs for shaping the radiation beam, and structure for generating a preliminary leaf sequence 
from a fluence map which supplies a desired radiation intensity profile, wherein the preliminary 
leaf sequence minimizes machine on-time and is generated without any leaf movement 
constraints. The system also includes structure for imposing at least one leaf movement 
constraint on the preliminary leaf sequence, and structure for applying at least one constraint 
elimination algorithm for adjusting the preliminary leaf sequence to minimize violations of the 
constraint while providing the fluence map and minimized on-time. The multi-leaf collimator 
can be a segmented multi-leaf collimator or a dynamic multi-leaf collimator. The system can 
include structure for reducing a tongue-and-groove underdose. The constraint can comprise a 
minimum separation distance between paired leaves. The structure for applying the constraint 
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can comprise structure for modifying at least one leaf pair in the preliminary leaf sequence to 
form a modified leaf sequence, wherein the structure for modifying comprises structure for 
identifying and adjusting positions of leaves in the preliminary leaf sequence which violate the 
minimum separation distance to provide at least the minimum separation distance; structure for 
modifying at least one leaf pair in the modified leaf sequence to produce a further modified leaf 
sequence, the further modified leaf sequence providing the desired intensity profile; structure for 
examining the further modified leaf sequence for violations of the minimum separation distance, 
and structure for iteratively adjusting the leaf sequence if at least one violation of the minimum 
separation distance is identified. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[00017] A better understanding of the present invention can be obtained when the following 
detailed description is considered in conjunction with the following drawings in which: 
[00018] FIG. 1 is a diagram of a radiation treatment device according an embodiment of the 
invention. 

[00019] FIG. 2 is a diagram of a multi-leaf collimator according to an embodiment of the 
invention. 

[00020] FIG. 3 is a flowchart illustrating a method according to an embodiment of the 
invention. 

[00021] FIG. 4 shows the geometry and coordinate system used in the Examples provided. 
[00022] FIG. 5(a) shows a piecewise continuous function while FIG. 5(b) shows the 
corresponding discretized intensity profile. 

[00023] FIG. 6 illustrates an exemplary leaf trajectory during SMLC delivery. 

[00024] FIG. 7(a) shows an intensity profile, while FIG. 7(b) shows the corresponding plan 

obtained using the algorithm SINGLE? AIR. 

[00025] FIG. 8 shows is leaf movement profile demonstrating bi-directional movement. 
[00026] FIG. 9 shows an MLC having three pairs of leaves. 

[00027] FIG. 10 illustrates the results of a leaf sequence modification process which 
eliminates minimum separation violations. 

[00028] FIGS. 11 (a) and (b) show intensity profiles to be delivered by adjacent leaf pairs t 
and / + 1 ^ respectively. 
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[00029] FIG. 12 shows the plans for/,(x) and from the profiles shown in FIGs. 1 1(a) 

and (b) obtained using algorithm MULTIPAIR. 

[00030] FIG. 13 illustrates a leaf trajectory plan obtained using Algorithm DMLC- 
SINGLEPAIR. 

[00031] FIG. 14 shows a profile using Algorithm DMLC-SINGLEPAIR. 
[00032] FIG. 15 shows a SMLC plan of a feasible profile. 

[00033] FIG. 16 shows the DMLC plan corresponding to the SMLC plan shown in FIG. 15 

using Algorithm DMLC-MINSINGLEPAIR. The plan is infeasible. 

[00034] FIG. 17 shows a bi-directional leaf movement profile. 

[00035] FIG. 18 illustrates "type 1" violations. 

[00036] FIG. 19 illustrates "type 2" violations. 

[00037] FIG. 20 is a cross section of two adjacent leaves. 

[00038] FIG. 21 shows a beams-eye view of the region to be treated by two adjacent leaf 
pairs. 

[00039] FIG. 22 shows a tongue-and-groove violation. 
[00040] FIG. 23 shows another tongue-and-groove violation. 
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DETAILED DESCRIPTION OF THE INVENTION 
[00041 1 A method of delivering radiation treatment to a patient using multi-leaf coUimation 
includes the steps of providing a radiation fluence map which includes an intensity profile. The 
fluence map is converted into a preliminary leaf sequence. The preliminary leaf sequence 
minimizes machine on-time and is generated without leaf movement constraints. The leaf 
movement constraint is imposed on the preliminary leaf sequence. At least one constraint 
elimination algorithm is then applied, the algorithm adjusting the preliminary leaf sequence to 
minimize violations of the constraint while providing the desired fluence map and a minimized 
radiation on-time. 

[00042] Turning now to FIG. 1, a block diagram of the radiation treatment device 100 
according to an embodiment of the invention is shown. An electron beam 105 is generated by an 
electron accelerator 106. The electron accelerator 106 includes an electron gun 1 10, a wave 
guide 112, and an evacuated envelope or guide magnet 1 13. A triggering system 1 14 generates 
injector trigger signals and supplies them to an injector 115. Based on the trigger signals, the 
injector 1 1 5 generates injector pulses which are fed to the electron gun 1 10 in the accelerator 106 
which results in the generation of electron beam 105. 

[00043] The electron beam 105 is accelerated and guided by wave guide 1 12. A high 
frequency signal source (not shown) is also provided, which supplies RF signals for the 
generation of an electromagnetic field which is supplied to wave guide 112. The electrons 
injected by the injector 115 and emitted by the electron gun 1 10 are accelerated by the 
electromagnetic field in the wave guide 112 and exit at the end opposite to electron gun 1 10 in 
electron beam 105. The electron beam 105 then enters guide magnet 113 and fi-om there is 
guided through window 1 1 7 along axis 118. After passing through a first scattering foil 1 19, the 
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beam goes through an opening 120 of a shield block 122 and encounters a flattening filter 123. 
Next, the beam is sent through a measuring chamber 125 in which the dose is determined. If the 
scattering foil 1 19 is replaced by a target, the radiation beam is an X-ray beam. Li this case, the 
flattening filter 123 may be absent. 

[00044] Beam shielding device is provided in the path of beam 105, comprising a plurality of 
opposing plates 131 and 132, only two of which are illustrated for convenience, hi one 
embodiment, other pairs of plates (not shown) are arranged perpendicular to plates 131 and 132. 
The plates 131 and 132 are moved with respect to axis 118 by a drive unit 134 to change the size 
and shape of the irradiated field. The drive unit 134 includes an electric motor which is coupled 
to the plates 131 and 132 and which is controlled by a motor controller 140. Position sensors 144 
and 145 are also coupled to the plates 131 and 132, respectively for sensing their positions. As 
noted above, the plate arrangement may alternatively include a multi-leaf collimator having a 
plurality of radiation blocking leaves. 

[00045] The motor controller 140 is coupled to a dosing unit 146 which includes a dosimetry 
controller and which is coupled to a central processing unit 148 for providing set values for the 
radiation beam for achieving given isodose curves. The output of the radiation beam is measured 
by a measuring chamber 125. In response to the deviation between the set values and the actual 
values, the dose control unit 146 supplies signals to a trigger system 114 which changes in a 
known manner the pulse repetition fi-equency so that the deviation between the set values and the 
actual values of the radiation beam output is minimized. In such a radiation device, the dose 
delivered is dependent upon movement of the collimator leaves 131 and 132. 
[00046] The central processing unit 148 is typically programmed by the therapist according to 
the instructions of an oncologist which performs beam optimization according to the present 
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invention so that the radiation treatment device carries out the prescribed radiation treatment 
while generally maximizing MU efficiency. The delivery of the radiation treatment is generally 
input through a keyboard 151, or other suitable data entry device. The central processing unit 
148 is further coupled to a dose control unit 146 that generates the desired values of radiation for 
controlling trigger system 1 14. The trigger system 1 14 then adapts the pulse radiation frequency 
and other parameters in a corresponding, conventional manner. The central processing unit 148 
further includes a control unit 156 which controls execution of the software and the opening and 
closing of the collimator plates 131 and 132 according to the present invention to deliver 
radiation according to a desired intensity profile generally having optimal MU efficiency. A 
monitor 160 is also provided. 

[00047] Figure 2 is a diagram of a multi-leaf collimator 200 according to an embodiment of 
the invention. Opposing leaf pairs 131a-n, 142a-n, each include a motor or drive unit 243a-n, and 
247a-n, respectively. The drive units are controlled by a control unit, such as control unit 156 
shown in Figure. 1 which controls execution of the software and the opening and closing of the 
collimator plates 131a-n and 142a-n according to the present invention to deliver radiation 
according to a desired intensity profile generally having optimal MU efficiency. 
[00048) Figure 3 is a flowchart illustrating a method 300 according to an embodiment of the 
invention. A radiation fluence map is first provided in step 301. In step 302 the fluence map is 
converted into a preliminary leaf sequence. The preliminary leaf sequence minimizes machine 
on-time and is generated without any leaf movement constraints. In step 303 at least one leaf 
movement constraint is imposed on the preliminary leaf sequence. In the case of a segmented 
multileaf system (SLMC) , the constraint includes a minimum leaf separation constraint, while in 
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the case of a dynamic multiieaf collimator system (DMLC), the constraint includes an 
interdigitation constraint. 

[00049] A tongue and groove constraint may also be applied to the preliminary leaf sequence, 
or the MU optimized sequence provided by method 300 without initially applying the tongue and 
groove constraint. A constraint elimination algorithm is then applied in step 304. This algorithm 
adjusts the preliminary leaf sequence to minimize violations of the constraint identified by the 
algorithm then further adjusts the leaf sequence to correct for errors induced in the intensity 
profile by corrections for the imposed constraint. In step 305 the algorithm checks whether the 
resulting leaf sequence provide violations below a predetermined limit, which can be no 
violations at all. If the resulting leaf sequence provides violations below a predetermined limit, 
the method is complete and a final leaf sequence is generated in step 306. If the resulting leaf 
sequence does not provides violations below a predetermined limit, step 304 is reapplied and the 
results rechecked in step 305, and so on iteratively, until no violations result. The resulting 
sequence provides the desired fluence map and minimizes radiation on-time. 
[00050] In the case of a segmented multiieaf system, an algorithm such as the SMLC 
MINSEPARATION algorithm described in the Examples operates by first solving the problem 
without any constraint to provide a preliminary leaf sequence that delivers the required intensity 
of radiation. This can be done using an algorithm such as the algorithm MULTIPAIR described 
in the Examples. This algorithm has been shown to theoretically generate an optimal MU 
efficiency leaf sequence without any constraint for SMLC systems. 

[00051] An algorithm then applies at least one leaf constraint. In this case, a minimum leaf 
pair separation constraint is generally applied. The algorithm MINSEPARATION can be used 
beginning with the preliminary leaf sequence determined above. It is possible that the 
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preliminary leaf sequence has violations of the minimum separation constraint, since the 
constraint was not considered when the preliminary leaf sequence was determined. Violations of 
the minimum separation constraint refer to time(s) in the sequence that opposing adjacent leaves 
in a given leaf pair are prescribed to be positioned too close together than is allowed. 
[00052] The minimum spacing depends on the particular model of the multileaf collimator 
used, and is generally specified by the manufacturer. Some models have a minimum separation 
constraint with minimum required spacing typically of the order of about 1 cm. Some models do 
not have a minimum separation constraint. 

[00053] Next, the positions of leaves in this sequence are scrutinized, such as starting fi-om 
left to right leaf pairs progressively. When a given position at which the minimum separation 
constraint is violated, the constraint elimination algorithm works to eliminate this violation. 
Corresponding leaf positions in the violating leaf pair(s) are adjusted to insure that the violating 
leaves are kept at least as far apart as is required to maintain minimum separation, while all other 
non- violating leaf positions are left unchanged. However, assuming some leaf positions are 
modified from the original preliminary leaf sequence, the new leaf sequence may no longer 
deliver the desired intensity profile as dictated by the fluence map provided. Accordingly, fiirther 
modifications to leaf pairs involved in the violation are generally required to maintain the same 
intensity profile as before. Once these modifications are applied, one constraint would be 
eliminated from violation from the original sequence, yet the same intensity profile will be 
maintained. 

[00054] Next, further violations of minimum separation in the modified leaf sequence that 
result from adjustments to recover the intensity profile as described above are identified and 
corrected. The process of correction for the intensity profile, then the leaf separation constraint, 
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and so on, are iteratively repeated until all violations of the minimum separation constraint are 
eliminated. 

[00055] It is possible that at some point in time, an intermediate leaf sequence generated is 
such that the leaves of a single leaf pair come too close together. In this case no feasible leaf 
sequence that delivers the exact desired intensity profile while also satisfying the minimum 
separation constraint may be possible. If this happens, the occurrence of which is preferably 
checked after every leaf modification to recover the intensity profile, the algorithm is terminated 
and no leaf sequence is generated. In such cases it is not possible to generate a leaf sequence that 
delivers the exact desired intensity profile. In this case, the inventive algorithm will not generate 
a leaf sequence as it will terminate without generating a profile. However, those having ordinary 
skill in the art will realize that a modified leaf sequencing algorithm can be designed to provide 
an approximate profile that can be generated in cases where an exact profile is not possible. 
[00056] If the intermediate leaf sequence generated is such that the leaves of a single leaf pair 
never come too close together, a leaf sequence is generated in which there are no further 
violations and which (theoretically) delivers the exact desired intensity. This sequence provides 
the optimal MU leaf sequence. 

[00057] The resulting leaf sequence has optimal MU efficiency because whenever a minimum 
separation constraint violation is detected, it is eliminated. The modifications applied to the leaf 
sequence to correct the intensity profile are such that they result in minimum loss of MU 
efficiency as possible for elimination of that particular instance of constraint violation. In other 
words, any other mechanism for elimination of that violation would result in higher loss of MU 
efficiency than the modification applied to recover the intensity profile. Since this is generally 
true every time the algorithm to correct for the intensity profile is applied, the final leaf sequence 



{WP151215;3| 



16 



has minimum loss of MU efficiency, thus providing a maximum MU efficiency or being optimal 

in MU efficiency. 

[00058] In the case of a dynamic multileaf system (DMLC), an algorithm, such as 
DMLC_INTERDIGITATION described in the Examples is applied without leaf constraints. 
This algorithm is the generally the same as the algorithm MINSEPARATION described for 
SMLC systems. However, the algorithm DMLC_MULTIPAIR is used which is a version of an 
algorithm MULTIPAIR applicable to dynamic multileaf coUimation (DMLC). The constraint 
used for DMLC systems is referred to as the interdigitation constraint (see Examples), as 
opposed to the minimum separation constraint described relative to the SMLC systems. The 
resulting modifications applied in this case are thus analogous, although different as compared to 
modifications for SMLC based systems according to the invention. 

[00059] For DMLC systems, there is no need to check whether a feasible leaf sequence results 
since the algorithm provided generally always results in a feasible optimal leaf sequence. So leaf 
sequence adjustments for the interdigitation constraint followed by adjustments to recover the 
intensity profile are applied iteratively until no violations remain. At this point the algorithm 
terminates and provides a leaf sequence with optimal MU efficiency. 
[00060] The invention can also be applied to correct for the tongue-and-groove effect. The 
tongue-and-groove effect results because in most multileaf collimators, the cross section of the 
leaves is not entirely rectangular. Rather, there is a tongue and groove overlap arrangement at the 
interface between adjacent leaves. The regions of the patient's body that fall below the tongue 
and groove overlap regions typically receive less radiation than is desired because these regions 
are shielded from the radiation when either of the two adjacent leaves shields the area adjacent to 
it. In some cases it is possible that some of these regions may not receive any radiation at all 
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though a certain amount of radiation has been prescribed to be delivered at these. To 
substantially reduce the resulting tongue-and-groove underdose, the invention provides 
algorithms to sequence the leaves as described below and in the Examples. 
[00061] The working of Algorithm TONGUEANDGROOVE described in the examples is 
generally the same as the algorithm MINSEPARATION disclosed herein for SMLC systems 
according to the invention, except for certain minor differences. In this case, the constraint being 
considered is the tongue-and-groove constraint, and not the minimum separation constraint as 
above to form the preliminary leaf sequence. The working of the algorithm 
TONGUEANDGROOVE-ID described in the Examples is the same as Algorithm 
MINSEPARATION except the constraint considered is the tongue-and-groove-id constraint 
(tongue-and-groove and interdigitation) and not the minimum separation constraint as above. 
The leaf sequence is then iteratively modified between adjustments for violations of the tongue- 
and-groove constraint and then violations of the resulting intensity profile analogous to the 
iterations described above for SMLC and DMLC systems above until no violations remain. A 
feasible leaf sequence generally always results. The resulting leaf sequence thus provides 
optimal MU efficiency. 



Examples 

[00062] It should be understood that the examples and embodiments described herein are 
for illustrative purposes only and that various modifications or changes in light thereof will be 
suggested to persons skilled in the art and are to be included within the spirit and purview of this 
application. The invention can take other specific forms without departing fi-om the spirit or 
essential attributes thereof 
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Example 1 

Leaf Sequencing Algorithms for Segmented Multileaf CoUimation 

[00063] The geometry and coordinate system used in this Example (and the other Examples) 
are shown in Figure 4. Delivery of profiles that are piecewise continuous are considered. Let 
I(x) be the desired intensity profile. The profiles first discretized to obtain the values at sample 
points Xq.x^^x^, I{x) is assigned the value I{x^) for < x < a:.^, , for each /.Now, /(x^) 

is the desired intensity profile. Figure 5(a) shows a piecewise continuous function while FIG. 
5(b) shows the corresponding discretized profile. The discretized profile can be efficiently 
delivered with the SMLC method. However, a SMLC sequence can be transformed to a dynamic 
leaf sequence by allowing both leaves to start at the same point and close together at the same 
point, so that they sweep across the same spatial interval. 

[00064] In this analysis it is assumed that the beam delivery begins when the pair of leaves is 
at the left most position. The initial position of the leaves is at Xq . Figure 6 illustrates the leaf 

trajectory during SMLC delivery. Let /,(jc.) and I.ix^) respectively denote the amount of 

Monitor Units (MUs) delivered when the left and right leaves leave position x. . Consider the 

motion of the left leaf. The left leaf begins at Xq and remains here until /^(x^) MUs have been 

delivered. At this time the left leaf is moved to , where it remains until //(x,) MUs have been 

delivered. The left leaf then moves to x^ where it remains until I^ix^) MUs have been delivered. 

At this time, the left leaf is moved to x^ , where it remains until //(x^) MUs have been delivered. 

The final movement of the left leaf is to x^ , where it remains until //(x^) = I^^ MUs have been 
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delivered. At this time the machine is turned off The total therapy time, TT{I^J^) , is the time 
needed to deliver /^^ MUs. The right leaf starts at ; moves to when I^{x^) MUs have 
been delivered; moves to x^ when /^(xj MUs have been delivered and so on. Note that the 
machine is off when a leaf is in motion. 

[00065] The following observations can be made. All MUs that are delivered along a 
radiation beam along ;c. before the left leaf passes x>^ fall on it. The greater the x value, the later 

the leaf passes that position. Therefore //(jc.) is a non-decreasing function. 

[00066J All MUs that are delivered along a radiation beam along x- before the right leaf 
passes X- , are blocked by the leaf The greater the x value, the later the leaf passes that position. 
Therefore I.^x^) is also a non-decreasing function. 

[00067] From these observations it is noted that the net amount of MUs delivered at a given 
point is given by Ii{x-)- I^{x-) , which must be the same as the desired profile I{x^) . 

[00068] When the movement of leaves is restricted to only one direction, both the left and 
right leaves move along the positive jc direction, from left to right (FIG. 4). Once the desired 
intensity profile, /(jc.) is known, one problem becomes that of determining the individual 

intensity profiles to be delivered by the left and right leaves, and such that: 

/(^/) = //(^/)-/.(^/XO</<m 

[00069] Where (//,/J is referred to as the treatment plan (or simply the plan) for / . Once 
the plan is obtained, it is possible to determine the movement of both left and right leaves during 
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the therapy. For each / , the left leaf can be allowed to pass X; when the source has delivered 
I,(xi) MUs. Also, the right leaf can be allowed to pass x, when the source has delivered 
MUs. In this manner the unidirectional leaf movement profiles for a plan are determined. 

[00070] From Equation 1, one way to determine /, and from the given target profile / is 
to begin with //(x^) = I(Xo) and /,(Xo) = 0 ; examine the remaining jc, s from left to right; 
increase /, whenever / increases; and increase whenever / decreases. Once /, and are 
determined the leaf movement profiles are obtained as explained in the previous section. The 
resulting algorithm SINGLEPAIR is shown below. Figure 7(a) shows the intensity profile while 
FIG. 7(b) shows the corresponding plan obtained using this algorithm. 
Algorithm SINGLEPAIR 

//(^o) = ^(^o) 
/.(^o) = 0 
For y = 1 to m do 
If (/(.:,)> /(X,.,) 
I,{Xj) = I,(x._,) + IiXj)-Iix..,) 

/,(x.) = /,(x,.,) 
Else 

I^Xj) = I^(Xj_,) + IiXj)-IiXj.,) 

End for 

[00071 ] Theorem 1 . It can be shown that the algorithm SINGLEPAIR obtains plans that are 
optimal in therapy time. 

[00072] Coroilary 1. Let /(x,) , 0 < i ^ m be a desired profile. Let /,(jc,) , and /,(x,) , 
0 < /■ < m be the left and right leaf profiles generated by Algorithm SINGLEPAIR. /,(x,.) and 
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/^(x.) , 0<i<m can be shown to define optimal therapy time unidirectional left and right leaf 
profiles for /(x.) , 0<i<m. The optimal therapy profiles proof follows from Theorem 1 . 

[00073) In the remainder of this example, (//,/,) is the optimal treatment plan for the desired 
profile / . 

[00074] The following observations are made about the optimal treatment plan (/p/^) 
generated using Algorithm SINGLEPAIR. 

[00075] Lemma 1. At each x- at most one of the profiles and changes (increases). 

[00076] Lemma 2. Let (Z^,/;?) be any treatment plan for / . 

(a) A(x,) = - I^ix.) = I,ix,)-IM) ^ 0,0 < / < m . 

(b) A(x.) is a non-decreasing function. 

[00077] Theorem 2. If the optimal plan (//,/^) violates the minimum separation constraint, 
then there is no plan for / that does not violate the minimum separation constraint. 

(000781 The separation between the leaves is determined by the difference in x values of the 
leaves when the source has delivered a certain amount of MUs. The minimum separation of the 
leaves is the minimum separation between the two profiles. This is referred to as the minimum 
separation . When the optimal plan obtained using Algorithm SINGLEPAIR is delivered, 

the minimum separation is . 

[00079] Corollary 2. Let S^j_^i„^^p,^ be the minimum leaf separafion in the plan (Ifyl,)^ Let 
^ud-min ^he minimum leaf separation in any (not necessarily optimal) given unidirectional 
plan. S^^_„,,„ < *S'„rf_;„/„(op/) • 
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[00080] Beam delivery when bi-directional movement of leaves is permitted is now examined. 
Whether relaxing the unidirectional movement constraint helps improve the efficiency of 
treatment plan is considered. 

[00081] For a given leaf (left or right) movement profile any x -coordinate can be classified as 
follows. Draw a vertical line at . If the line cuts the leaf profile exactly once x is called an 
unidirectional point of that leaf profile. If the line cuts the profile more than once, x is a bi- 
directional point of that profile. A leaf movement profile that has at least one bi-directional point 
is a bi-directional profile. All profiles that are not bi-directional are unidirectional profiles. Any 
profile can be partitioned into segments such that each segment is a unidirectional profile. When 

the number of such partitions is minimal, each partition is called a stage of the original profile. 

J- 

Thus unidirectional profiles consist of exactly one stage while bi-directional profiles always have 
more than one stage. 

[00082] Now referring to FIG. 8, the leaf movement profile, , shows the position of the left 
leaf as a function of the amount of MUs delivered by the source. The leaf starts from the left 
edge and moves in both directions during the therapy. Clearly, is bi-directional. The 

movement profile of this leaf consists of stages S^,S^ and . In stages 5, and the leaf 
moves fi-om left to right while in stage the leaf moves firom right to left. x. is a bi-directional 
point of B^ . The amount of MUs delivered at x^ is I, + . In stage 5, , when amount of 
MUs have been delivered, the leaf passes x- . Now, no MU is delivered at Xj till the leaf passes 
over X. in . Further, MUs are delivered to x^ in stages and 53 . Thus, Ii{Xj) = L^ +^3. 
Here, L^^I^.L^-I^-I^^ a:^ is a unidirectional point of B^ . The MUs delivered at x,^ are 
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L^-I^> Note that the intensity profile /, is different from the leaf movement profile , unlike 
in the unidirectional case. 

[000831 Lemma 3. Let (//,/,) be a plan delivered by the bi-directional leaf movement 
profile pair {B^,B^) (i.e., B^ and B^ are, respectively, the left and right leaf movement profiles) 

(a) /, is non-decreasing. 

(b) is non-decreasing. 

[00084] From Lenmia 3 it is noted that a bi-directional leaf movement profile B delivers a 
non-decreasing intensity profile. This non-decreasing intensity profile can also be delivered 
using a unidirectional leaf movement profile as described earlier. This profile is referred to as the 
unidirectional leaf movement profile that corresponds to the bi-directional profile B and will be 
denoted by U to emphasize that it is unidirectional. Thus, every bi-directional leaf movement 
profile has a corresponding unidirectional leaf profile that delivers the same amount of MUs at 
each X. as does the bi-directional profile. 

[00085] Theorem 3. The unidirectional treatment plan constructed by Algorithm 
SINGLEPAIR is optimal in therapy time even when bi-directional leaf movement is permitted. 

[00086] Let and be unidirectional leaf movement profiles that deliver the desired 
profile /(jc.) . Let B^ and B^ be a set of bi-directional left and right leaf profiles such that and 

correspond to B^ and B^ respectively, i.e., (5^, 5J delivers the same plan as (C/pi7J.The 
minimum separation of leaves in this bi-directional plan is referred to herein as {B^.B^) . 

[00087] Theorem 4. i'^^.,,,/^ < S^^_^^^ for a bi-directional leaf movement profile pair and its 
cori'esponding unidirectional profile. 
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[00088] Theorem 5. If the optimal unidirectional plan violates the minimum 

separation constraint, then there is no bi-directional movement plan that does not violate the 
minimum separation constraint. 

[00089] In this example portion, a single pair of leaves is used to deliver intensity profiles 
defined along the axis of the pair of leaves. However, in a real application, it is required to 
deliver intensity profiles defined over a 2-D region. Multi-Leaf Collimators (MLCs) are used to 
deliver such profiles. An MLC is composed of multiple pairs of leaves with parallel axes. 
Figure 9 shows an MLC that includes three pairs of leaves - (Zl,/?l),(Z,2,i22) and (Z3,/23). 
I1,L2,I3 are left leaves and Rl,R2,R3 are right leaves. Each pair of leaves is controlled ' 
independently. If there are no constraints on the leaf movements, the desired profile is divided 
into a set of parallel profiles defined along the axes of the leaf pairs. Each leaf pair / then 
delivers the plan for the corresponding intensity profile /.(x) . The set of plans of all leaf pairs 
forms the solution set. This set is referred to as the treatment schedule (or simply, the schedule). 

[00090] In practical situations, however, there are some constraints on the movement of the 
leaves. As shown above, the minimum separation constraint requires that opposing pairs of 
leaves be separated by at least some distance ( 5^,.^ ) at all times during beam delivery. In MLCs 

this constraint is applied not only to opposing pairs of leaves, but also to opposing leaves of 
neighboring pairs. For example, in FIG. 9, LI and iil, L2 and /?2, L3 and R3, LI and R2, 
LI and R\ , L2 and R3 , L3 and R2 are pairwise subject to the constraint. The term intra-pair 
minimum separation constraint is used to refer to the constraint imposed on an opposing pair of 
leaves and inter-pair minimum separation constraint to refer to the constraint imposed on 
opposing leaves of neighboring pairs. Recall that, as shown above, it was proved that for a single 
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pair of leaves, if the optimal plan does not satisfy the minimum separation constraint, then no 
plan satisfies the constraint. In this section an algorithm is presented to generate the optimal 
schedule for the desired profile defined over a 2-D region. The algorithm is then modified to 

generate schedules that satisfy the inter-pair minimum separation constraint. 

[00091] An optimal schedule without the minimum separation constraint is now described. 
Assume n pairs of leaves. For each pair, there are m sample points. The input is represented as 
a matrix with n rows and m columns, where the / th row represents the desired intensity profile 
to be delivered by the / th pair of leaves. The Algorithm SINGLEPAIR is applied to determine 
the optimal plan for each of the n leaf pairs. This method of generating schedules is described in 
Algorithm MULTIPAIR shown below. 
Algorithm MULTIPAIR 
For(/ = l;/<w;/ + +) 

Apply Algorithm SINGLEPAIR to the / th pair of leaves to obtain plan *at delivers the 

intensity profile I^{x) . 
End For 

[00092] Lemma 4. Algorithm MULTIPAIR generates schedules that are optimal in therapy 
time. 

[00093] An optimal algorithm with inter-pair minimum separation constraint is now 
described. The schedule generated by Algorithm MULTIPAIR may violate both the intra- and 
inter-pair minimum separation constraints. If the schedule has no violations of these constraints, 
the schedule is the desired optimal schedule. If there is a violation of the intra-pair constraint, 
then it follows from Theorem 2 that there is no schedule that is fi-ee of constraint violation. So, 
assume that only the inter-pair constraint is violated. All violations of the inter-pair constraint 
are eliminated starting from the left end, i.e., from . To eliminate the violations, those plans of 
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the schedule that cause the violations are modified. The schedule is scanned from along the 
positive X direction looking for the least x^ at which is positioned a right leaf (say Ru) that 
violates the inter-pair separation constraint. After rectifying the violation at x^ with respect to 
Ru, other violations are looked for. Since the process of eliminating a violation at x^ , may at 
times, lead to new violations at Xj.Xj <x^/\t is sometimes needed to retract a certain distance 
(shown as 5^.^ ) to the left, every time a modification is made to the schedule. The scanning and 
modification process is then restarted from the new position. The process continues until no 
inter-pair violations exist. Algorithm MINSEPARATION below outlines this procedure. 
Algorithm MINSEPARATION 
Assume no intra-pair violations exist 

i. x = Xq 

ii. While (there is an inter-pair violation) do 

iii. Find the least x^, x^>x, such that a right leaf is positioned at x^ and this right leaf has an 
inter-pair separation violation with one or both of its neighboring left leaves. Let u be the 
least integer such that the right leaf Ru is positioned at x^ and Ru has an inter-pair 
separation violation. Let Lt denote the left leaf (or one of the left leaves) with which Ru has 
an inter-pair violation. Note that / € {m - 1, w + 1} . 

iv. Modify the schedule to eliminate the violation between Ru and Lt . 

V. If there is now an intra-pair separation violation between Rt and Lt , no feasible schedule 
exists, terminate. 

vi. a: = jc -iS^. 

V nttn 

vii. End While 
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[00094] Let M = {{luJxMhnhX-^ihnJnr)) be the schedule generated by Algorithm 
MULTIPAIR for the desired intensity profile. 

Let A^(p) = ((/w,,/..p),(/2/,.A.pX'-(^/p.^«.p)) be the schedule obtained after step (iv) of the 
Algorithm MINSEPARATION is applied p times to the input schedule M. Note that 

M = A^(0) . 

To illustrate the modification process to eliminate violations, FIG. 10 is provided. To 
simplify, only two neighboring pairs of leaves are shown in FIG. 10. Suppose that the (/? + i)th 
violation occurs when the right leaf of pair u is positioned at and the left leaf of pair 
^ / e {w - 1, w + 1} , arrives at , x^ - x^ < 5^, . Let x^=x^- . To remove this inter-pair 
separation violation, (IffpJtrp) jg modified ^^^^^ profiles of N(p) are not modified. The 
new I^^p (i.e., /,/(^+]) ) is defined below. 

max{I,/^(x),/,,(x) + A/} x^<x<x^ 
where ^ = LpM^ ^A^u) = h^^i- Ar(p.i)W = WoW^^rW^ where I,{x) is the target 
profile to be delivered by the leaf pair / . Since //^(^+,) differs fi'om /^^ for x>x^=x^- S^-^ there 
is a possibility that N(p-^l) has inter-pair separation violations for right leaf positions 
x>x^=x^- 5^.^ . Since none of the other right leaf profiles are changed firom those of N(p) 
and since the change in /^^ only delays the rightward movement of the left leaf of pair t , no 
inter-pair violations are possible in A^(/7 + l) for x<x^ ^^v^^min- One may also verify that 
since and /,^o are non-decreasing functions of ;c , so also are /,,^ and , p>0. 
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[00095) Lemma 5. Let F = {{luJxMhnliX-^ilnnhr)) be any feasible schedule for the 
desired profile, i.e., a schedule that does not violate the intra- or inter-pair minimum separation 
constraints. Let S{p) , be the following assertions. 
ll,{x)>l,^{x), 0<i<n,x,<x<x^ 
ii,/- W>/.^(x), 0<i<n,x,<x<x^ 
S{p) is true for p>0, 

[00096] Lemma 6. If an intra-pair minimum separation violation is detected in step (v) of 
MINSEPARATION, then there is no feasible schedule for the desired profile. 

[00097] Suppose that there is a feasible schedule F and that leaf pair t has an intra-pair 

minimum separation violation in N(p\p > 0 . From Lemma 5 it follov^s that 

(a) r^i(x)>I^i^{xlx^<x<x^„ 

(b) I^^(x)>I,^{xlx,<x<x„, 

[00098] An instance is now described where an inter-pair minimum separation violation is 
detected in step (v) of the algorithm MINSEPARATION. Figure 1 1(a) and (b) show intensity 
profiles to be delivered by adjacent leaf pairs t and / + 1 ^ respectively. The plans for/^(x) and 
If^^ix) are obtained using algorithm MULTIPAIR are shown in Figure 12. Each of these plans 
{(I^i{x)J^^{x)) and /(,+,);. W)) is feasible, i.e., there is no intra-pair minimum 

separation ( = 7 ). However, when MINSEPARATION is applied (for simplicity consider 
leaf pairs t and / + 1 in isolation), it detects an inter-pair minimum separation violation between 
and 7,^ , when 7^,^,), arrives at ;c = 6 and 7^ is positioned at ;c = 1 1 . To eliminate this 
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violation, is positioned at jc = 4 (since 1 1 - 4 = 7 = S'^j, ) and its profile is raised from 
x = 4 . Consequently is also raised from ^ = 4 resulting in the plan W»-'(,+i)ri 
This modification results in an intra-pair violation for pair / + 1 , when Z^,^,)^, is at x = \ and 
is at X = 4 . From Lemma 6, there is no feasible schedule. 

[00099] For N{p\p > 0 and every leaf pair j,\<j<n, define 
^Jip(^-^) = ^..(^-i) = 0, A = -/,^(x,.,),0 < / < m and 

^jrp(^i) = ^rp(^i) " ^rp(^i-i)'0 - ^ ' • Notice that Ajipix^) gives the time (in monitor units) for 
which the left leaf of pair j stops at position jc,. . Let Aj^^ix.) and Aj^^ix.) be zero for all x- 
when y = 0 as well as when y = /i + 1 . 

[000100] Lemma 7. For every j\l< j<n and every i,\<i<m, 

)} (2) 



[000101] Lemma 8. For every y,l< j<n and every i,l</<m, 

A,,^(;c,) = A,,^(^,)-(/,(x,)-//x,.,) (3) 

where /^.(jc_,) = 0. 

[000102] Notice that once a right leaf u moves past a:^ , no separation violation with respect to 
this leaf is possible. Therefore, (see algorithm MINSEPARATION) <a:^. Hence, 
A < A .;o(Jc,) , and A < A ..^(x,),^, - < x, < < y < « . Starting with these upper 
bounds, which are independent of /? , on Aj^p{x.) , x^ - iS^.^ < x,. < x^ and using Equations 2 and 
3, an upper bound on the remaining A^.,^(;C;)s and e/^a^^(;f,.)s are computed (from right to left). 
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The remaining upper bounds are also independent of p . Let the computed upper bound on 

^jipi^i) be Uji{x.) . It follows that the therapy time for (IjipJjrp) is at most 
'^maxU) = Xo^,'5m ^A^i^ ' Thcrcforc, the therapy time for N(p) is at most 

[000103] Theorem 6. The following are true of Algorithm MINSEPARATION: 

a. The algorithm terminates. 

b. When the algorithm terminates in Step (v), there is no feasible schedule. 

c. Otherwise, the schedule generated is feasible and is optimal in therapy time. 

[000104] Corollary 3. When S^^„ = 0 , the Algorithm Minseparation always generates an 
optimal feasible schedule. 

Example 2 

Optimal Sequencing of Dynamic Multileaf Collimators 

[000105] The geometry and coordinate system used in this Example are the same as shown in 
FIG. 4. The delivery of intensity map produced by the optimizer is considered. It is important to 
note that the intensity map from the optimizer is always a discrete matrix. The spatial resolution 
of this matrix is similar to the smallest beamlet size. The beamlet size typically ranges from 5-10 
mm. Most leaf sequencing algorithms for DMLC delivery published in literature and 
implemented in commercial software first convert the discrete intensity matrix into a piecewise 
continuous function with linear interpolation between the intensity values thus yielding smooth 
and continuous intensity map. This is believed to be unnecessary if the leaf sequencing algorithm 
can exactly reproduce the optimized intensity map. The leaf sequencing algorithm described 
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herein is designed to reproduce exactly the same intensity map as it comes out of the optimizer. 
This ensures an exact reproduction of the planned dose distribution. With this consideration, let 

^(•^^ be the desired intensity profile along ^ axis. The discretized profile from the optimizer 

XX X 

gives the intensity values at sample points . . . ' " . 

[000106] It is assumed that the sample points are uniformly spaced and that 

Ax = X,,, - X,, 0 < / < m /(X) i5 ^ggigj^^d the value ^^^^ > for " ^ < , for each i . Now, ^^^^ > 

is our desired intensity profile, i.e., ^^'^'^ is a measure of the number of MUs for which ^' , 

0 < / < m ^ needs to be exposed. 

[000107] Figure 5(b) shows a discretized profile, which is the output from the optimizer at 

discrete sample points ^o^-^i^ . . . ' . In this Example, this profile is delivered with Dynamic 
Multileaf CoUimation (DMLC). 

[000108] In this analysis it will be assumed that ^^^o) > 0 ^j^^ ^(^«) > ^ and that when the 
beam delivery begins the leaves can be positioned anywhere. It is also assumed that the leaves 

can move with any velocity ^ , "^'"^ " ^ " ^""^"^ , where ^'"^ is the maximum allowable velocity 
of the leaves. Figure 13 illustrates an exemplary leaf trajectory during DMLC delivery. In this 

I (x) I (x) 

Example, the leaves move from left to right. Let and , respectively, denote the 

amount of Monitor Units (MUs) delivered when the left and right leaves leave position ^' . 

X I (x ) 

Consider the motion of the left leaf The left leaf begins at ° and remains here until MUs 

have been delivered. At this time the left leaf leaves and is moved to ,where it remains until 

1 (x) X I (x ^ 

'"^ MUs have been delivered. The left leaf then moves to ^ where it remains until 
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X I (x ) 

MUs have been delivered. At this time, the left leaf is moved to ^ ^ where it remains until 

X I (x ^ 

MUs have been delivered. Then it moves to ^ , where it remains until «^ MUs have been 

delivered. The final movement of the left leaf is to . The left leaf arrives at when 
MUs have been delivered. At this time the machine is turned off The total therapy time, 

TTU^J^) ^ ji^g ^jj^^^ needed to deliver ^'^"^'^ MUs. The right leaf starts at and begins to move 

rightaway till it reaches ; leaves when ^'•^^2) mus have been delivered; leaves when 

r\ a) y[\js have been delivered, and so on. Note that the machine is on throughout the 

X 

treatment. All MUs that are delivered along a radiation beam along ' before the left leaf passes 
^' fall on it. Similarly, all MUs that are delivered along a radiation beam along ^' before the 
right leaf passes ^' , are blocked by the leaf So the net amount of MUs delivered at a point is 
given by ^/(-^i) "^'•(^z) , which must be the same as the desired profile ^^"^'^ . 

[000109] Theorem 1. The following are true of every leaf pair trajectory that delivers a 

discrete profile: <^ 

(a) The left leaf must reach at some time. 

(b) The right leaf must reach at some time. 

(c) The left leaf must reach at some time. 

(d) The right leaf must reach at some time. 

[0001 10] As noted earlier, the velocity of leaves cannot exceed some maximum limit (say ) 
in practice. This implies that the leaf profile cannot be horizontal at any point. From FIG. 13, it 
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can be seen that the time needed for a leaf to move from ^' to is " ^i )/^/mu , if O is the 
flux density of MUs from the source, the number of MUs delivered in this time along a beam is 
^ . So, ^ <l>*(^,>i -^,)/v.a. = Phi^^^Ka. . The same is true for 

the right leaf profile . 

[0001 11] An optimal unidirectional algorithm for one pair of leaves is now described. When 
the movement of leaves is restricted to only one direction, both the left and right leaves move 

I(x ) • 

along the positive ^ direction, from left to right. Once the desired intensity profile, ^ is 
known, the problem becomes that of determining the individual intensity profiles to be delivered 

by the left and right leaves, and such that: 

/(x,) = //W-/.(^/XO</</w 

[000112] and are subject to the maximum velocity constraint. ^^^'^^^ is referred to as the 
treatment plan (or simply plan) for ^ . Once the plan is obtained, the movement of both left and 
right leaves during the therapy will be determined. For each ^ , the left leaf can be allowed to 

pass ' when the source has delivered MUs. Also, it is possible to allow the right leaf to 

X I (x) 

pass ' when the source has delivered MUs. In this manner, unidirectional leaf movement 

profiles for a plan are determined. Some observations about the leaf profiles are made below. 

[000113] Theorem 2. In every unidirectional plan the leaves begin at and end at . A 
proof of this theorem follows from Theorem 1 and the imidirectional constraint. 
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X 

[0001 14] Lemma 1. In the region between each pair of successive sample points, say ' and 



, both leaf profiles maintain the same shape, i.e., one is merely a vertical translation of the 
other. 

[0001 151 Lemma 2. In an optimal plan, both leaves must move at between every 
successive pair of sample points they move across. 

[000116] Corollary L In an optimal plan, no leaf stops at an ^ that is not one of the ^' s. 

From Equation 1, it can be seen that one way to determine and from the given target 

profile ^ is to begin from set ^/(^o) = ^(^o) and ^'•(•^o) = 0. examine the remaining ^' s to 

the right; increase ' at ' whenever ^ increases and by the same amount (in addition to the 

minimum increase imposed by the maximum velocity constraint); and similarly increase 

whenever ^ decreases. This can be done till is reached. So the treatment begins with the 
leaves positioned at the leftmost sample point and ends with the leaves positioned at the 

rightmost sample point. Once and are determined the leaf movement profiles are obtained 

as explained earlier. Note that leaves are moved at the maximum velocity ^""^ whenever they are 
to be moved. The resulting Algorithm DMLC-SINGLEPAIR is shown below. Figure 14 shows a 

profile / and FIG. 13 shows the corresponding plan ^^t'^r) obtained using Algorithm DMLC- 
SINGLEPAIR. 

Algorithm DMLC-SINGLEPAIR 

/.(^o) = 0 
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For J to f" do 

(/(;c,)>/(^,.,)) 

/, (Xj ) = /, (V. ) + li^j ) - /(V, ) + ^ * ^/v_ 
Else 

/,(x.) = /,(x^..,) + <D*Ax/v„„ 
End for 

[0001 1 7] It has been shown that the Algorithm DMLC-SINGLEPAIR obtains plans that are 
optimal in therapy time. However, the proof is complex and indirect. A simpler and direct proof 
is provided below. 

(0001 1 81 Theorem 3. Algorithm DMLC-SINGLEPAIR obtains plans that are optimal in 
therapy time. 

[000119] Corollary 2. Let ^^^'^ , 0 ^ / ^ w be a desired profile. Let ^'^^'^ and ^'^^'^ , 

0 < /■ < ra be the left and right leaf profiles generated by Algorithm DMLC-SINGLEPAIR. ^'^^'^ 

and ^'■^''^'^ , 0<i<m define optimal therapy time unidirectional left and right leaf profiles for 

^i^i) 0<i<m 

[000120] In the remainder of this section, ^^i'^--^ is the optimal treatment plan generated by 

Algorithm DMLC-SINGLEPAIR for the desired profile ^ . 
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[000121] The following observations are made about the optimal treatment plan ^ '^^ 
generated using Algorithm DMLC-SINGLEPAIR. 

[000122] Lemma 3. At most one of the leaves stops at each ^' . 

[000123] Lemma 4. Let ^^^'^^^ be any treatment plan for I . 

A(^,) = I,{x,)-I,{x,) = /^(x,) - > 0,0 < / < /w 

(b) ^^"^'^ is a non-decreasing function. 

[000124] Corollary 3. Let ^'^^'^ (^r(^i)) and ^L^^i^ ^^^f^^'^^ respectively, denote the 
amount of time for which the left (right) leaf stops at ^' in plans ^^i'^r) and (^^'^^^ . Then 

(a) A.(^,)^A/U,), 

(b) ^r{x,)>\{x.) 

[000125] Some MLCs have a minimum separation constraint that requires the left and right 

e 

leaves to maintain a minimum separation at all times during the treatment. Notice that in the 
plan generated by Algorithm DMLC-SINGLEPAIR, the two leaves start and end at the same 

point. Accordingly, they are in contact at and . When ^, the minimum separation 

constraint is violated at and "^^ . In order to overcome this problem, the Algorithm DMLC- 
SINGLEPAER is modified to guarantee minimum separation between the leaves in the vicinity of 

the end points ("^^ and ). In particular, the left leaf is permitted to be initially positioned at a 

point -^0 "^min and the right leaf to be finally positioned at ^ '^^min ^ jhe only changes 

made relative to Algorithm DMLC-SINGLEPAIR are for the movement of the left leaf from 
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to ^0 and for the right leaf from ^™ to ^'"^ . The movement of the left leaf from to (and a 

symmetric definition for the right leaf from to ^'"') is defined to be such that it maintains a 

distance of exactly from the right leaf at all times. Once the left leaf reaches ^ it follows 
the trajectory as before. While this modification results in the exact profile being delivered 

between ^ and it also results in some undesirable exposure to the intervals ^ o'* o^' and 

i^my^m') In the remainder of this section exposures of this kind permissible will be considered 

X X 

permissible, provided the exact profile is delivered between o and . Note that for most 

(x X J (x X * } 

commercially available accelerators, undesirable exposure to the intervals ^ and ^ '"'^ 

X X 

can be avoided by positioning the back-up jaws at ^ and respectively. However, a difficulty 
arises when the number of monitor units delivered at the time the left leaf reaches in this new 

(TT\ I (x ^ XX 

plan (call it ^ '^^) is greater than . This would prevent using the old plan from ^ to « , 

since the leaf cannot pass "^^^ before it arrives there. However, that if the left leaf were to arrive at 

^0 any earlier, it would be too close to the right leaf. In the discussion that follows it is shown 
that in this and other cases where the original plan violates the constraint, there are no feasible 

solutions that deliver exactly the desired profile between and , while permitting exposure 
outside this region. The modified algorithm, which is referred to herein as DMLC- 
MESfSINGLEPAIR, is described below. Note that the therapy time of the plan produced by 
DMLC-MINSINGLEPAIR is the same as that of the plan produced by DMLC-SINGLEPAIR. 
Therefore, the modified plan is optimal in therapy time. 
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Algorithm DMLC-MINSINGLEPAIR 

1 . Apply Algorithm DMLC-SINGLEPAIR 

2. Modify the profile of the left leaf from to and the right leaf from ^« to to maintain 
a minimum interleaf distance of ^""^^ . Call this profile ^^^'^r) , 

X I (x ) 

3. If the number of MUs delivered when the left leaf arrives at «> is greater than there is 
no feasible solution. End. 

4. Else output '^'^ There is a feasible solution only if ^^^'^r) feasible. 

[000126] Theorem 4. (a) ""^^^max > or (b) If the plan ^^^'^'•^ generated by DMLC- 
MINSINGLEPAIR violates the minimum separation constraint, then there is no plan for ^ that 
does not violate the minimum separation constraint. 

[000127] The separation between the leaves is determined by the difference in ^ values of the 
leaves when the source has delivered a certain amount of MUs. The minimum separation of the 
leaves is the minimum separation between the two profiles. This minimum separation is referred 

c 

to herein as . When the optimal plan obtained using Algorithm DMLC-SINGLEPAIR is 

S 

delivered, the minimum separation is . 

[000128] Corollary 4. Let ^"d-mimopt) the minimum leaf separation in the plan ^^''^'-^ . Let 

c 

"^ud-min be the minimum leaf separation in any (not necessarily optimal) given unidirectional 
c < c 

[000129] Theorem 5. If Algorithm DMLC-MDSFSINGLEPAIR terminates in Step 3, then there 
is no plan for ^ that does not violate the minimum separation constraint. 
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The DMLC leaf sequence is now compared with small sequencing. 

[000130] Lemma 5. Let the minimum separation between the leaves in the optimal SMLC plan 

be . Let the minimum leaf separation in the plan generated by algorithm DMLC- 

MINSINGLEPAIR be . Then " . 

The following result immediately follows and can be easily verified. 

(000131) Corollary 5, All profiles that have feasible plans using DMLC have feasible plans 
using SMLC. There exist profiles for which there are feasible plans using SMLC, but no feasible 
plans using DMLC. 

[000132] Figure 15 shows a SMLC plan for a feasible profile. The corresponding DMLC plan 
obtained using Algorithm DMLC-MINSINGLEPAIR shown in Figure 16 is seen to be 
infeasible. 

[000133] Beam delivery when bi-directional movement of leaves is permitted using the 
invention is now described. The effects of relaxing the unidirectional movement constraint on the 
efficiency of treatment plan is considered. For a given leaf (left or right) movement profile any 
^ -coordinate is classified as follows. Draw a vertical line at ^ . If the line cuts the leaf profile 
exactly ^ is referred to as a unidirectional point of that leaf profile. If the line cuts the profile 
more than once, ^ is referred to as a bi-directional point of that profile. A leaf movement profile 
that has at least one bi-directional point is a bi-directional profile. All profiles that are not bi- 
directional are unidirectional profiles. Any profile can be partitioned into segments such that 
each segment is a unidirectional profile. When the number of such partitions is minimal, each 
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partition is called a stage of the original profile. Thus unidirectional profiles consist of exactly 
one stage while bi-directional profiles always have more than one stage. 

D 

[000134] Figure 1 7 shows a bi-directional leaf movement profile, ^ , which shows the 
position of the left leaf as a function of the amount of MUs delivered by the source. The 

s s s s s 

movement profile of this leaf consists of stages ' ' ^ and 3 ^ stages > and ^ the leaf 

moves from left to right while in stage 2 the leaf moves from right to left. ^ is a bi-directional 

B I B 
point of ^ . Let ^ be the intensity profile corresponding to the leaf movement profile ^ . 

^i(^) gives the number of MUs delivered at ^ using movement profile . The amount of MUs 

delivered at ^ is ' 2 stage ' , when » amount of MUs have been delivered, the leaf 

X X X S I — / 

passes ' . Now, no MU is delivered at ^ till the leaf passes over ^ in 2 . Further, ^ 2 jynjs 

are delivered to ^> in stages and . Thus, ^' ) = A + ^2 ^ ^^ere " and ^2 = A ^ A . 

X B X L = I 

* is a unidirectional point of ^ . The MUs delivered at * are ^ 4 Note that the intensity 

r D 

profile ^ is different from the leaf movement profile ^ , unlike in the unidirectional case. 

[000135] Lemma 6. Let and be the intensity profiles corresponding to the bi-directional 
leaf movement profile pair ^ ^' (i.e., ^ and ^ are, respectively, the left and right leaf 
movement profiles). Let ^^^'^ ^ ^"^'^ ~ ^^'^ , 0 < ^ ^ w , be the intensity profile delivered by 
(^^'^^^Then 

(b) ^/r(^/.l)^^/?(-^/) + ^*^/^mav. 
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(B B ) 

[000136] From Lemma 6 it is noted that every bi-directional leaf movement profile ^ 

delivers an intensity profile ^^^'^R^ that satisfies the maximum velocity constraint. Hence, 

Uu^r) is deliverable using a unidirectional leaf movement profile. This profile is referred to 
herein as the unidirectional leaf movement profile that corresponds to the bi-directional profile. 
Thus every bi-directional leaf movement profile has a corresponding unidirectional leaf profile 

that delivers the same amount of MUs at each as does the bi-directional profile. 

[000137] Theorem 6. The unidirectional treatment plan constructed by Algorithm DMLC- 
SINGLEPAIR is optimal in therapy time even when bi-directional leaf movement is permitted. 

[000138] Incorporation of the minimum separation constraint will now be described. Let ^' 

II I(x ) B 

and ' be unidirectional leaf movement profiles that deliver the desired profile ^ . Let ' 

and be a set of bi-directional left and right leaf profiles such that and correspond to ^' 

and respectively, i.e./^''^^^ delivers the same plan as ^^I'^r) ^ This is called the minimum 

separation ofleaves in this bi-directional plan ^ "^bd-mm ^ ^^ud-mm is the minimum separation 

of leaves in ^^i'^r) 

[000139] Theorem 7, ^rf-""" ~ "^ud-min for every bi-directional leaf movement profile pair 
i^n^r) and its corresponding unidirectional profile ^i'^r) ^ 

[000140] Theorem 8, If the optimal unidirectional plan ^^^'^r) violates the minimum 
separation constraint, then there is no bi-directional movement plan that does not violate the 
minimum separation constraint. 
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[000141] In a practical application, it is needed to deliver intensity profiles defined over a 2-D 
region. Multi-Leaf Collimators (MLCs) are used to deliver such profiles. An MLC is composed 
of multiple pairs of leaves with parallel axes. Referring again to FIG. 9, FIG. 9 shows an MLC 
that has three pairs of leaves - (A^l),(^2,/?2) (I3,/?3) L\,L2,L3 are left leaves and 
R\,R2,R3 ^^^^ leaves. Each pair of leaves is controlled independently. If there are no 
constraints on the leaf movements, The desired profile is divided into a set of parallel profiles 
defined along the axes of the leaf pairs. Each leaf pair ^ then dehvers the plan for the 

corresponding intensity profile \ The set of plans of all leaf pairs forms the solution set. 
This set is referred to as the treatment schedule (or simply schedule). 

[000142] In practical situations, however, there are constraints on the movement of the leaves. 
As discussed earlier, the minimum separation constraint requires that opposing pairs of leaves be 

separated by at least some distance ( ) at all times during beam delivery. In some MLCs this 
constraint is applied not only to opposing pairs of leaves, but also to opposing leaves of 
neighboring pairs. For example, in FIG. 9, ^1 and ^1 , ^2 and ^2 ^ L3 and ^3 ^ LI and ^2 ^ 
L2 and ^1, ^2 and and ^2 are pairwise subject to the constraint. The term intra-pair 

minimum separation constraint is used to refer to the constraint imposed on an opposing pair of 
leaves and inter-pair minimum separation constraint to refer to the constraint imposed on 
opposing leaves of neighboring pairs. The inter-pair minimum separation constraint with 
S =0 

is of special interest and is referred to as the interdigitation constraint. Recall that, 
above it was proven that for a single pair of leaves, if the optimal plan does not satisfy the 
minimum separation constraint, then no plan satisfies the constraint. In this section an algorithm 
is disclosed which generates the optimal schedule for the desired profile defined over a 2-D 
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region. The algorithm can be modified to generate schedules that satisfy the interdigitation 
constraint. It is noted that in the discussion on single pair of leaves, it was assumed that 

^(•^o) > 0 and that ^^^"'^ ^ ^ . However, with multiple leaf pairs, the first and last sample points 
with non-zero intensity levels could be different for different pairs. Hence, this assumption will 
no longer be made.. 

[000143] Assume ^ pairs of leaves. For each pair, there are ^ sample points. The input is 
represented as a matrix with ^ rows and ^ columns, where the ' th row represents the desired 
intensity profile to be delivered by the ^ th pair of leaves. The Algorithm DMLC-SENGLEPAIR 
is applied to determine the optimal plan for each of the ^ leaf pairs. This method of generating 

schedules is described in Algorithm DMLC-MULTIPAIR shown below. Note that since 

X X X ' X * 

are not necessarily non-zero for any row, ^ is replaced by ' and is replaced by * m 

Algorithm DMLC-SINGLEPAIR for each row, where and ^« , respectively, denote the first 
and last non-zero sample points of that row. Also, for rows that contain only zeroes, the plan 

simply places the corresponding leaves at the rightmost point in the field (call it ). 

(000144] Lemma 7. Algorithm DMLC-MULTIPAIR generates schedules that are optimal in 
therapy time. 

Algorithm DMLC-MULTIPAIR 
Pqj.(/ = !;/<«;/ + +) 

Apply Algorithm DMLC-SINGLEPAIR to the ' th pair of leaves to obtain plan ^^'''^^^ that 

/ (x) 

delivers the intensity profile ^. 
End For 
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[000145] The schedule generated by Algorithm DMLC-MULTIPAIR may violate the 

S =0 

interdigitation constraint. Note that no intra-pair constraint violations can occur for 
Accordingly, the interdigitation constraint is essentially an inter-pair constraint. If the schedule 
has no interdigitation constraint violations, it is the desired optimal schedule. If there are 
violations in the schedule, all violations of the interdigitation constraint can be eliminated 

starting from the left end, i.e., from . To eliminate the violations, plans of the schedule that 

cause the violations are modified. The schedule from o along the positive ^ direction is 

scarmed looking for the least at which is positioned a right leaf (say ^« ) that violates the 

X R 

inter-pair separation constraint. After rectifying the violation at with respect to « other 

violations are looked for. Since the process of eliminating a violation at ^ , may at times, lead to 

new violations involving right leaves positioned at , it is preferred to search afresh from "^^ 
every time a modification is made to the schedule. The scanning and modification process is 
continued until no interdigitation violations exist. Algorithm DMLC-INTERDIGITATION 
shown below discloses the procedure. 
Algorithm DMLC-INTERDIGITATION 

i. -^"^0 

ii. While (there is an interdigitation violation) do 

iii. Find the least ^ such that a right leaf is positioned at and this right leaf has an 
interdigitation violation with one or both of its neighboring left leaves. Let " be the least 

D X R 

integer such that the right leaf " is positioned at ^ and " has an interdigitation violation. 
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T D 

Let ' denote the left leaf with which " has an interdigitation violation. Note that 
r € {w + 1} case has violations with two adjacent left leaves, let ^ = " . 



R L 

iv. Modify the schedule to eliminate the violation between " and ' . 



V. »' 



vi. End While 

[0001461 Let ^ ^ ((^i/'^ir)»(A/,/2rX --^C^./' A,r)) thc schedule generated by Algorithm 
DMLC-MULTIP AIR for the desired intensity profile. 

[000147] Let ^^^^ " ((^up^^irpWiipyhrpl - yi^nip^hrp)) schedule obtained after Step (iv) 
of Algorithm DMLC-INTERDIGITATION is applied P times to the input schedule ^ . Note 
that^ = ^(0). 

[000148] To illustrate the modification process other examples are presented. There are two 
types of violations that may occur. They are called Type 1 and Type 2 violations and the 
corresponding modifications Type 1 and Type 2 modifications herein. To simplify, only two 

neighboring pairs of leaves are shown. Suppose that the th violation occurs between the 

right leaf of pair " , which is positioned at , and the left leaf of pair ^'^ ^ {w -l,w + 1} ^ 

[000149] In a Type 1 violation, the left leaf of pair ^ starts its sweep at a point ^^^^^^(^>P) ^ 
(see FIG. 18). To remove this interdigitation violation, modify ^^^'p'^'^/'^ to (^'/(p+O'^/rCP+i)) as 
follows. The leaves of pair ^ are allowed to start at and move them at the maximiun velocity 
towards the right, till they reach ^Start{t,p) ^j^^ number of MUs delivered when they 



V 
max 
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reach xStart(t,p) A . Raise the profiles ^""^''^ and ^^^'^'«''^(^/'), by an amount 

/, = <i>*{xStart{t,p)-xyv^ ^j^^ ^^^^j^. 

= x,)/v^ ^,<x< xStart{t, p) 

l,ipix)-^h X >xStart(t,p) 

I/r(p+i)W = A/(p+i)W-^W^^here ^'^^^ is the target profile to be delivered by the leaf pair K 

[000150] A Type 2 violation occurs when the left leaf of pair ^ , which starts its sweep from 

^ - passes before the right leaf of pair " passes (Figure 19). In this case, '^"(''*'> is as 
defined below. 

where ^-^un>M~^iip^^v) = h~h Qnce again, ^'r(P*\)^^^~ ^ii(p+i)(x)~^t(x) ^ where '^'^''^^ is 
the target profile to be delivered by the leaf pair ^ . 

[000151] In both Type 1 and Type 2 modifications, the other profiles of ^^P^ are not 
modified. Since ^"(p*') differs from ^"p for ^"^^ there is a possibility that ^^P"^^^ has inter- 
pair separation violations for right leaf positions ^ - ^ Since none of the other right leaf 
profiles are changed from those of ^^P^ and since the change in only delays the rightward 
movement of the left leaf of pair ^, no interdigitation violations are possible in ^^P'^^^ for 
^ "^^ It can also be verified that that since ^''^ and ^""^ are feasible plans that satisfy the 
maximum velocity constrains, so are and ^^'^ , P>0 
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[000152] 



Lemmas. I jn>^^StartU.p)) = 0 ^ i<j<m^ p>0. 



[000153] Corollary 6. A Type 2 violation in which ^"p^^"^ ~ ^ cannot occur. 
[000154] Lemma 9. In case of a Type 1 violation, ^^"p'^'^p^ is the same as ^^no'^iro) . 

[000155] Note that ^""^^^ and ^"^^^^ are defined only for ^ - ^-^^^^'(^ . In the sequel, the 

convention that ^ ~ ^ 4p(^)) jg true is adopted whenever ^""^^^ (^/rp W ) is undefined, 

irrespective of whether ^ is defined or not. 

[000156] Lemma 10. Let ^ ~ (('^i"Ar)»(A/>Ar)» - '(^B/'^r)) be any feasible unidirectional 
schedule for the desired profile, i.e., a unidirectional schedule that does not violate the 
interdigitation constraint. Let ^^P^ , be the following assertions. 

a. W ^ litM) 0<i<n,x,<x<x„ 

b. 0<i<n,x,<x<x„, 

^^P^ is true for /'^O. 

[000157] Lemma 11. For the execution of Algorithm DMLC-INTERDIGITATION 

(a) ^^"^^ Type 1 violations can occur. 

(b) Type 2 violations can occur. 

T 

(c) Let be the optimal unidirectional therapy time for the input matrix. The time complexity 
0{mn + n*min{nm,T^}) 
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[000158] Note that Lemma 1 1 provides two upper bounds of on the complexity of Algorithm 
DMLC-INTERDIGITATION: and ^^^*'^^^'"'^'"«'^^ . In most practical situations, 

^max < and so * ^^^^'^'^-aJ) can be considered a tighter bound. 

[0001 591 Theorem 9. The following are true of the Algorithm DMLC-INTERDIGITATION: 

a. The algorithm terminates. 

b. The schedule generated is feasible and optimal in therapy time unidirectional schedule. 



Example 3: 

Optimal Leaf Sequencing for Elimination of Tongue-and-Groove Underdosage 

[000160] The delivery of optimal radiation intensity map produced by the inverse treatment 
planning system is considered in this Example. The intensity matrix from the optimizer generally 
has a spatial resolution that is similar to the smallest beamlet size. The beamlet size generally 

ranges from 5-10 mm. ^^^^ is the desired intensity profile. The discretized profile from the 

X X X X 

optimizer gives the intensity values at sample points °' ^ ' . . . ' . It is assumed that the 

sample points are uniformly spaced and that ^ ~ - x.,0<i<m I(x) jg assigned the value 

^(^/) for ^/ ^ ^ < ^i^i ^ for each ^ , Now, ^^^'^ is the desired intensity profile, i.e., ^^^'^ is a 

measure of the number of MUs for which , 0<i<m^ needs to be exposed. Figure 5(b) is 
referred to again which shows a discretized profile, which is the output from the optimizer. This 
discretized profile is delivered either with the Segmental Multileaf CoUimation (SMLC) method 
or with Dynamic Multileaf CoUimation (DMLC). An SMLC sequence can be transformed to a 
dynamic leaf sequence by allowing both leaves to start at the same point and close together at the 
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same point, so that they sweep across the same spatial interval. The current Example is 
developed for SMLC delivery. 

[000161] In this analysis it is assumed that the beam delivery begins when the pair of leaves is 

at the left most position. The initial position of the leaves is in this paper, it is assumed that 
leaves may move only from left to right. That is unidirectional leaf movement is assumed. 
Figure 6 is again referred to which illustrates the leaf trajectory during SMLC delivery. Let 

^'^■^'^ and ^'^^'^ respectively denote the amount of Monitor Units (MUs) delivered when the 

left and right leaves leave position ^' . Consider the motion of the left leaf. The left leaf begins at 

and remains here until ^^^^o) ]v[Us have been delivered. At this time the left leaf is moved to 

, where it remains until ^'^^i^ MUs have been delivered. The left leaf then moves to where 

it remains until ^/^-^s) mus have been delivered. At this time, the left leaf is moved to , where 

it remains until ^i^^^^ MUs have been delivered. The final movement of the left leaf is to ^'^ , 

where it remains until ^^(-^7) - Liax jsAJJs have been delivered. At this time the machine is turned 

off. The total therapy time, ^^^^I'^r) ^ the time needed to deliver ^'"^^ MUs. Note that the term 

therapy time is used to refer to the beam-on time. The right leaf starts at 2 ; moves to ^ when 

/ / y \ X I (X ) 

'^^2/ have been delivered; moves to ^ when ''^ MUs have been delivered and so on. 
Note that the machine is off when a leaf is in motion. The following observations are made: 

X X 

i. All MUs that are dehvered along a radiation beam along ' before the left leaf passes ' 

I (x) 

fall on it. Greater the ^ value, later the leaf passes that position. Therefore is a 

non-decreasing function. 
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ii. All MUs that are delivered along a radiation beam along ^' before the right leaf passes 

, are blocked by the leaf. Greater the ^ value, later the leaf passes that position. 

J (x) 

Therefore is also a non-decreasing function. 

[000162] From these observations, it can be seen that the net amount of MUs delivered at a 
point is given by ^'^"^'^ " ^ which must be the same as the desired profile ^^^'^ . 

[000163] Once the desired intensity profile, ^^^'^ is known, the single leaf pair problem 
becomes that of determining the individual intensity profiles to be delivered by the left and right 

leaves, and such that: 

/W = //W-/.(^/),0</<m ... 



[000164] (^"^'^^ is referred to as the treatment plan (or simply plan) for ^ . Once the plan is 

obtained, it is possible to determine the movement of both left and right leaves during the 

X I (x) 

therapy. For each ^ , the left leaf can be allowed to pass ' when the source has delivered 

X I (x) 

MUs. Also, the right leaf can be allowed to pass ' when the source has delivered MUs. 

A plan is unidirectional if and are unique, ^ " i.e., each leaf passes over 

each point only once. The Algorithm SINGLEPAIR is presented below which an be used to 
obtain a unidirectional plan. 

[000165] Algorithm SINGLEPAIR 

/,K) = A^o) 
/.(^o) = 0 
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For to f" do 

^^(I{Xj)>I{XjJ 

Else 

fMj) = IMj.,) + IiXj.,)-IiXj) 
I,{Xj) = I,iXj_,) 

End for 

[000166] Theorem L Algorithm SINGLEPAIR obtains plans that are optimal in therapy time 
even when bidirectional leaf movement is permitted. 

1000167] Corollary 1. Let ^^^'^ , 0 ^ ^ be a desired profile. Let ^'^^'^ and ^'•^^'^ , 

0 < / < AW be the left and right leaf profiles generated by Algorithm SINGLEPAIR. ^'^^'^ and 

^r(^i) ^ 0<i<m define optimal therapy time unidirectional left and right leaf profiles for ^^^'^ , 
0</<7 

[000168] Lemma 1. Let ^^^'^^^ be any treatment plan for I . 

(a) ^(^/) = hi^i) - hi^i) = h i^i) ' tMi) > 0,0 < /• < //i 

(b) ^^^'^ is a non-decreasing function. 

[000169] In a real applications, it is necessary to deliver intensity profiles defined over a 2-D 
region. Multi-Leaf Collimators (MLCs) are used to deliver such profiles. An MLC is composed 
of multiple pairs of leaves with parallel axes. Figure 9 is again referred to which shows an MLC 
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that has three pairs of leaves - (^l^l),(^2,/?2) (13, /?3) I1,Z,2,Z3 j^f^ j^^^gg 

R\,R2,R3 j^gj^^ leaves. Each pair of leaves is controlled independently. If there are no 
constraints on the leaf movements, desired profile is divided into a set of parallel profiles defined 
along the axes of the leaf pairs. Each leaf pair ^ then delivers the plan for the corresponding 

intensity profile ^'^'^^ . The set of plans of all leaf pairs forms the solution set, referred to herein 
as the treatment schedule (or simply schedule). A schedule in w^hich all plans are unidirectional 
is a unidirectional schedule. Only unidirectional schedules are considered in this Example. 

[000170] Assume ^ pairs of leaves. For each pair, there are ^ sample points. The input is 
represented as a matrix v^ith ^ rov^s and ^ columns, where the ' th row represents the desired 
intensity profile to be delivered by the ' th pair of leaves. The Algorithm SINGLEPAIR is 
applied to determine the optimal plan for each of the ^ leaf pairs. This method of generating 
schedules is described in Algorithm MULTIPAIR shown below. Since the individual plans of the 
leaf pairs are optimal in therapy time, it can be shown the resulting schedule also is optimum in 
therapy time. 
Algorithm MULTIPAIR 
Pqj.(/ = 1;/<a2;/ + +) 

Apply Algorithm SINGLEPAIR to the ^ th pair of leaves to obtain plan ^^n'^^r) that dehvers the 

intensity profile ^ . 
End For 

Optimal algorithms with interdigitation and tongue-and-groove constraints are now presented. 
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[000171] In practical situations, there are some constraints regarding movement of the leaves. 
The minimum separation constraint requires that opposing pairs of leaves be separated by at least 

some distance ( ) at all times during beam delivery. In some MLCs this constraint is applied 
not only to opposing pairs of leaves, but also to opposing leaves of neighboring pairs. For 
example, in FIG. 9, L\ and ^1 , ^2 and ^2 , L3 and ^3 ^ LI and R2 , L2 and ^1 , L2 and 
R3 ^ L3 and ^2 are pairwise subject to the constraint. The term intra-pair minimum separation 
constraint is used herein to refer to the constraint imposed on an opposing pair of leaves and 
inter-pair minimum separation constraint to refer to the constraint imposed on opposing leaves of 

5 =0 

neighboring pairs. The inter-pair minimum separation constraint with is of special 

interest and is referred to herein as the interdigitation constraint. 

[000172] In most commercially available MLCs, there is a tongue-and-groove arrangement at 
the interface between adjacent leaves. A cross section of two adjacent leaves is .depicted in 
Figure 20. The width of the tongue-and-groove region is ^ . The area under this region gets 
underdosed due to the mechanical arrangement. Figure 21 shows a beams-eye view of the region 

to be treated by two adjacent leaf pairs, ^ and ^ + ^ . Consider the shaded rectangular areas ^^"^'^ 

and that require exactly ^'^^'^ and MUs to be delivered, respectively. The 

tongue-and-groove overlap area between the two leaf pairs over the sample point ' , ^»'+^^ *^ , is 

A (x) I (x) 

colored black. Let the amount of MUs delivered in be '-'^^^ . The following lemma 

(Lemma 2) is a consequence of the fact that is exposed only when both and 

/+iv ij are exposed. 
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[0001731 Leminal. 0</ <m, ^ 

[0001 74] Schedules in which ^' '^^ ^ ^ ™^ ^' ^'^^ are said to be free of tongue-and- 
groove underdosage effects are now described. Unless treatment schedules are carefully 

designed, it is possible that nun{/,(x.),/,^,(x.)} ^^^^ ' and ^ For example, in a 

schedule in which ^^rM = ^0^ W(^J = 50 V,,(^J = 60^ 

(^vv) = A/ (^.v) ' (^^) = 50 - 50 = 0 ^^^^ ^j^^^ ^j^g ^^g^^ 

min{/,(^J,/,,.(xJ} = /,,,(^J-/,(xJ = 60-50 = 10 ^^^^ ^^^.^ ^^^p^^ ^j^^^ 

could be 0 even when is arbitrarily large. 

[000175] The leaf sequencing algorithms described in Examples 1 (SMLC) and 2 (DMLC) 
generate schedules that satisfy the inter-pair minimum separation constraint and are optimal in 
therapy time. However, these algorithms do not account for the tongue-and-groove effect. In this 
Example, two algorithms are presented. Algorithm TONGUEANDGROOVE generates 
minimum therapy time unidirectional schedules that are free of tongue-and-groove underdosage. 
Algorithm TONGUEANDGROOVE-ID generates minimum therapy time unidirectional 
schedules that are free of tongue-and-groove underdosage while simultaneously satisfying the 
interdigitation constraint. 

[000176] The following lemma provides a necessary and sufficient condition for a 
unidirectional schedule to be free of tongue-and-groove underdosage effects. 

[000177] Lemma 3. A unidirectional schedule is free of tongue-and-groove underdosage 
effects if and only if. 
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For 0<i<m^ l<t <n 

[000178] Lemma 3 is equivalent to a statement that the time period for which a pair of leaves 

A (x) 

(say pair ^ ) exposes the region ^ is completely contained by the time period for which 

pair ^ 1 exposes region "^' '^^ ^"^'^ , or vice versa, whenever ^'^^'^ ^ and ^ ^ . Note that 

if either ^'^^'^ or is zero the containment is not necessary. This necessary and sufficient 

condition of Lemma 3 is referred to herein as the tongue-and-groove constraint condition. 
Schedules that satisfy this condition will be said to satisfy the tongue-and-groove constraint. 

[000179] Note that the schedule generated by Algorithm MULTIPAIR described above may 
violate the tongue-and-groove constraint. If the schedule has no tongue-and-groove constraint 
violations, it is the desired optimal schedule. If there are violations in the schedule, all violations 
of the tongue-and-groove constraint can be theoretically eliminated starting from the left end, 

i.e., from . To eliminate the violations, plans of the schedule that cause the violations are 

modified. The schedule is scanned from along the positive ^ direction looking for the least 

at which there exist leaf pairs " , ^ t e {u-\,u-\'\} ^ ^j^^j violate the constraint at . After 

rectifying the violation at , other violations are checked for. Since the process of eliminating a 

violation at ^*%may at times, lead to new violations at it is sometimes necessary to search 
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afresh from every time a modification is made to the schedule. However, it can be proven 

that a bound of ^^^^ on the number of violations that can occur at After eliminating all 

violations at a particular sample point, the analysis moves to the next point, i.e., increment 
^ and look for possible violations at the new point. The scanning and modification process 
continues until no tongue-and-groove constraint violations exist. Algorithm 
TONGUEANDGROOVE shown below outlines the procedure. 
Algorithm TONGUEANDGROOVE 

i ^ = ^0 

ii. While (there is a tongue-and-groove violation) do 

iii. Find the least ^^-^ ^ such that there exist leaf pairs " , " + 1 , that violate the tongue- 
and-groove constraint at 

iv. Modify the schedule to eliminate the violation between leaf pairs " and ^ + 1 . 

V. = 

vi. End While 

[000180] Let ^ ^ ((A/' Ar)'(^2/>^2r)» ">(^«/>0) be the schedule generated by Algorithm 

MULTIPAIR for the desired intensity profile. 

\ 

[000181] Let ^^^^ ^ iihip^hrp)^ihipJ2rp\-^ihipJnrp)) the schedule obtained after Step (iv) 
of the Algorithm TONGUEANDGROOVE is applied P times to the input schedule M . Note 
that^ = ^(0). 
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[000182] To illustrate the modification process examples are used. To simplify, only two 
neighboring pairs of leaves are shown. Suppose that the th violation occurs between the 

leaves of pair " and pair ^^^-^^ dX^^, Note that ^'^^ ^^"'^ ^"''^ ^^^^ , as otherwise, either (b) or 

(c) of Lemma 3 is true. In case ^^'^^^^^^ ^ hipM ^ g^^p w and ^ Now, ^''^^^"^ ^ hipM ^ ^^iQ 
sequel, these are referred to as " and ^ values as the " and ^ of Algorithm 
TONGUEANDGROOVE. From Lemma 3 and the fact that a violation has occurred, it follows 

that ^^'p^^'""^ ^ Lrp(^w) remove this tongue-and-groove constraint violation, ^^f^p'^^'p\s 
modified. The other profiles of ^^P^ are not modified. 

[000183] The new plan for pair ^ , (^'/(p+D'^/^/'+o) is as defined below. 

Kip (^iv ) ~ Kip i^w ) - Krp i^w ) ~ Krp i-^w ) jj^gn 

M = {lip (^) ^o^x<x^ 
^.(x) + A/ <a:<x 

where ^ = KipM- KipM ^/r(p+i)W = ^//(/;+i)W~^rW ^ ^here ^'^^^ is the target profile to be 

delivered by the leaf pair ^ . 

Otherwise, 

T/.„W + Ar' <x<x 

where ^ =hrpM-I,rpM -^//(p+dW = Ar(p+i)W + ^/W^ ^l^grg AW js the target profile to 
be delivered by the leaf pair ' . 
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[000184] . The former case is illustrated in Figure 22 and the latter case is illustrated in 
Figure 23. 

[000185] Since ^^'/(p+»'Ar(p+i)) differs from ^^^^p'^^^p) ^j. ^^^^^ is a possibility that 

^{p + 1) is involved in tongue-and-groove violations for ^ ^ . Since none of the other leaf 
profiles are changed from those of ^^P^ no tongue-and-groove constraint violations are possible 
in ^(^"*"^) for One may also verify that since ^''^ and are non-decreasing functions 

of ^,so also are ^^^^ and P^^, 

[000186] Lemma 4. Let ^ ^ ^^hiyhMhnhX' X^nn^nr)) be any unidirectional schedule for 

the desired profile that satisfies the tongue-and-groove constraint. Let '^^^^ be the following 
assertions. 

^ I,{x)>I,^{x)^ Q<i<n,x,<x<x^ 
4 W ^ 4p W ^0<i<n,x,<x< x^ 
^(P^ is true for /^^O 

Algorithms for the elimination of tongue-and-groove effect and interdigitation will now 
be described. As pointed out above, the elimination of tongue-and-groove constraint violations 
does not guarantee elimination of interdigitation constraint violations. Therefore the schedule 
generated by Algorithm TONGUEANDGROOVE may not be fi-ee of interdigitation violations. 
The algorithm described herein for obtaining schedules that simultaneously satisfy both 
constraints, Algorithm TONGUEANDGROOVE-ID, is similar to Algorithm 
TONGUEANDGROOVE. The only difference between the two algorithms lies in the definition 
of the constraint condition. To be precise, the following definition is provided: 
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Deflnition 1. A unidirectional schedule is said to satisfy the tongue-and-groove-id constraint if 
. /Jx,) < ^ ^ , or 

For 0<i<m^ \ <t <n 

[000187] The only difference between this constraint and the tongue-and-groove constraint is 
that this constraint enforces condition (a) or (b) above to be true at all sample points ' including 
those at which ^'^^'^ " ^ and/or " ^ . 

[000188] Lemma 5. A schedule satisfies the tongue-and-groove-id constraint if it satisfies the 
tongue-and-groove constraint and the interdigitation constraint. 

[000189] Algorithm TONGUEANDGROOVE-ID finds violations of the tongue-and-groove-id 
constraint fi-om left to right in exactly the same maimer in which Algorithm 
TONGUEANDGROOVE detects tongue-and-groove violations. Also, the violations are 
eliminated as before, i.e., as prescribed by Equations 2 and 3 and illustrated in FIGS. 22 and 23, 
respectively. Algorithm TONGUEANDGROOVE-ID is shown below. All notation used in the 
algorithm and the related discussion in the remainder of this example is also the same as that 
used above and corresponds directly to the usage in Algorithm TONGUEANDGROOVE. 

Algorithm TONGUEANDGROOVE-ID 

i. ^"^0 

ii. While (there is a tongue-and-groove-id violation) do 
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iii. Find the least ^^^^^-^ ^ such that there exist leaf pairs " , " ^ that violate the tongue- 

and-groove-id constraint at 

iv. Modify the schedule to eliminate the violation between leaf pairs " and " ^ . 

V. 

vi. End While 

[000190] In the remainder of this example the term "algorithm" is used to refer to the 
Algorithm TONGUEANDGROOVE or Algorithm TONGUEANDGROOVE-ID and "violation' 
is used to refer to tongue-and-groove constraint violation or tongue-and-groove-id constraint 
violation (depending on which algorithm is considered), unless explicitly mentioned. 

[000191] The execution of the algorithm starts with schedule ^ at ^ and sweeps to the 
right, eliminating violations from the schedule along the way. The modifications applied to 

eliminate a violation at prescribed by Equations 2 and 3, modify one of the violating profiles 

for ^ ~ ■^H' From the unidirectional nature of the sweep of the algorithm, it is clear that the 

modification of the profile for ^ ^ can have no consequence on violations that may occur at 

the point Therefore it suffices to modify the profile only at at the time the violation at 
is detected. The modification can be propagated to the right as the algorithm sweeps. This can be 

done by using an matrix ^ that keeps track of the amount by which the profiles have 

been raised. ^(■^ '^) denotes the cumulative amount by which the J th leaf pair profiles have 

been raised at sample point ^* from the schedule A/ generated using Algorithm MULTIPAIR. 

When the algorithm has eliminated all violations at each it moves to to look for possible 
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violations. It first sets the ^^'^'^^th column of the modification matrix equal to the ^th column 
to reflect rightward propagation of the modifications. It then looks for and eliminates violations 



at and so on. 



[000192] The process of detecting the violations at is described further. It is shown below 
that if one carefully selects the order in which violations are detected and eliminated, the number 

of violations at each 0 < w < w be ^^^^ . 

[000193] Lemma 6. The algorithm can be implemented such that ^^^^ violations occur at 
each 0<w<w. 

[000194] The bound is achieved using a two pass scheme at In pass one adjacent leaf pairs 
(l,2),(2,3),...,(w -1,«) ^j.^ checked, in that order, for possible violations at In pass two, 
violations are checked for in the reverse order, i.e., ('^-Uw),(«-2,w-l),...,(l,2) ^^^^ 
adjacent pairs ^^'^ l<i<n checked exactly twice for possible violations. It can be seen 
that if a violation is detected in pass one, either the profile of leaf pair ' or that of leaf pair ' "^^ 
may be modified (raised) to eliminate the violation. However, in pass two only the profile of pair 
' may be modified. This is because the profile of pair ' is not modified between the two times it 
is checked for violations with pair ^ ^ . The profile of pair ' ^ , on the other hand, could have 
been modified between these times as a result of violations with pair ' + 2 Therefore in pass 
two, only ' can be a candidate for ^ (where ^ is as explained in the algorithm) when pairs 
{hi + 1) are examined. From this it also follows that when pairs are subsequently 

examined in pass two, the profile of pair ^ will not be modified. Since there is no violation 
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between adjacent pairs (1' 2), (2, 3),.. .,(/,/ + 1) ^Y^^^ ^^^^ ^^^^ ^j^^g^ p^jj-g jg g^gj. 
examined again, it follows that at the end of pass two there can be no violations between pairs 
(^^ + 1) l<i<n 

[000195] Lemma 7. For the execution of the algorithm, the time complexity is . 
Follows from Lemma 6 and the fact that there are sample points. 

[000196] Theorem 2. 

a. Algorithms TONGUEANDGROOOVE and TONGUEANDGROOVE-ID terminate. 

b. The schedule generated by Algorithm TONGUEANDGROOVE is free of tongue-and- 
groove constraint violations and is optimal in therapy time for unidirectional schedules. 

c. The schedule generated by Algorithm TONGUEANDGROOVE-ID is free of 
interdigitation and tongue-and-groove constraint violations and is optimal in therapy time 
.for unidirectional schedules. 

[000197] It is to be understood that while the invention has been described in conjunction with 
the preferred specific embodiments thereof, that the foregoing description as well as the 
examples which follow are intended to illustrate and not limit the scope of the invention. Other 
aspects, advantages and modifications within the scope of the invention will be apparent to those 
skilled in the art to which the invention pertains. 
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